Phases of granular segregation in a binary mixture 
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We present results from an extensive experimental investigation into granular segregation of a 
shallow binary mixture in which particles are driven by frictional interactions with the surface of 
a vibrating horizontal tray. Three distinct phases of the mixture are established viz; binary gas 
(unsegregated), segregation liquid and segregation crystal. Their ranges of existence are mapped 
out as a function of the system's primary control parameters using a number of measures based 
on Voronoi tessellation. We study the associated transitions and show that segregation can be 
suppressed is the total filling fraction of the granular layer, C, is decreased below a critical value, 
Cc, or if the dimensionless acceleration of the driving, 7, is increased above a value 7c. 



I. INTRODUCTION 

Granular materials, i.e. ensembles of macroscopic dis- 
crete particles, are ubiquitous in our every day life, na- 
ture and of crucial importance in industrial processes 
P, |2- The study of granular media has a long tradi- 
tion amongst engineers and geologists who have had suc- 
cesses in specific problems using a combination of prac- 
tical experience and empirical knowledge. More recently 
the physics community has taken an increased interest 
in granular materials since they pose a number of funda- 
mental questions which challenge current ideas in non- 
equilibrium statistical mechanics Interactions be- 
tween granular particles are intrinsically dissipative since 
energy is lost due to both inelastic collisions and frictional 
contacts. Hence, any dynamical study of a granular en- 
semble requires an energy input which typically takes 
the form of vibration or shear In this sense, de- 
tailed investigations of granular media provide examples 
of canonical systems where dynamical processes are far 
from equilibrium |^. 

An interesting and counter-intuitive feature of particu- 
late matter is segregation of binary assemblies, where an 
initially uniform mixture of particles can spontaneously 
de-mix into its constituent components under flow Q- 
Typically, the species of particles may differ in size, den- 
sity, rigidity or surface properties. Such differences can 
often lead to separation and hence clustering of like par- 
ticles 0, 0. Intriguingly, segregation does not always 
happen and the conditions for its occurrence are difficult 
to predict. An extensive account of the issues involved 
can be found in the following reviews 0, Q • 

The phenomena has been recognized for a considerable 
period of time 0| but, despite more than half a century 
of research, the underlying nature of the mechanisms in- 
volved are not yet understood. Apart from posing vari- 
ous theoretical questions, insight into segregation would 



be beneficial for many applications. These include areas 
such as agriculture, geophysics, material science, and sev- 
eral branches of engineering, e.g. involving preparation 
of food, drugs, detergents, cosmetics, and ceramics 1^. 
In many of these examples processing and transporting of 




FIG. 1: The experimental set up. (a) Three dimensional 
schematic diagram of the experimental apparatus, (b) Typ- 
ical segregation patterns of the binary granular mixtures in 
the segregation liquid phase. Poppy seeds appear as dark 
grey regions and phosphor-bronze spheres as light gray. The 
frame was taken after 3min of vibration from an initially 
homogeneous mixture, (C = 0.708, F — 2, f = 12Hz, 
A = ±1.74mm). (c) Photographs of a phosphor-bronze 
sphere and (d) a poppy seed. 
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mixtures can lead to undesired separation. Segregation 
of mixtures can be brought about by the simple actions 
of pouring, shaking, vibration, shear and fluidization and 
is also found in industrial processes where the objective 
is to achieve particle mixing . 

Over the last decade, significant interest in segregation 
has arisen in the physics community. A number of small 
scale laborato ry expe riments have been reported on ver- 
tically [H HI, HIS and horizontally |23]llljllj2| vi- 
brated beds, filling and emptying of vessels ^4 [l2 5 |^ l26ll2 ^ 
and rotating cylindrical drums ,25. 28, 29, S^jISJ- The 
spatial distribution and dynamics of segregation of large 
and small particles appears to depend on a number of 
factors besides size difference including density ratio, fric- 
tion between particles, shape of boundaries, particle ve- 
locity and the effects of the interstitial air. A great deal 
of research is required on the topic since the parameter 
space of such a system is large and seemingly trivial de- 
tails turn out to have important effects in the segregation 
of mixtures. Designing simple and well controlled lab- 
oratory based experiments and developing appropriate 
models is therefore essential to make progress in under- 
stand segregation which in turn may give insights into 
the industrial problems outlined above. 

We have recently developed an experiment to study the 
segregation of shallow layers of granular binary mixtures 
under horizontal vibration [s^, Issl 133 . l35l |. The exis- 
tence and self-organisation of three phases of segregation 
by systematically starting from homogeneous binary mix- 
tures have been uncovered. These are, binary gas (unseg- 
regated), segregation liquid and segregation crystal and 
they exist over ranges of the total filling fraction of the 
layer, C. The principal result is the discovery of critical 
phenomena in the segregation process. This implies the 
existence of a transition point in C below which the layer 
remains mixed and above which segregation occurs. An 
overview of this work will be given in Section [llll 

In this paper we present the results of an extensive ex- 
perimental investigation where we show that this phase 
behavior of the granular mixture is robust over a range 
of control parameter space. Moreover, we investigate 
the role of the driving on the segregation process and 
uncover an additional transition between the segregated 
and mixed phases as the dimensionless acceleration of 
the tray is increased. We perform particle tracking and 
use the distribution of positions of the centers of one of 
the particle types to define a number of measures derived 
from Voronoi tessellation to characterize the state of the 
mixture. This detailed study points to a robust behavior 
of our granular mixture which invites the development of 
a predictive theoretical model. 

This paper is organized as follows. A description of the 
experimental apparatus is given in Section ^] In Section 
mil we describe the nature of segregation patterns ob- 
served in our system and briefiy review some results from 
our previous work. The microscopic measures obtained 
from Voronoi tessellation, namely the local Voronoi den- 
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cussed in Sections II V Al and II V Bl In SectiorTvlwe report 
the results of an exploration of the parameter space of 
the system including the aspect ratio of the cell (Section 
IV A|) , the filling fraction of the mixture fScction lV Bp and 
forcing parameters (Section IV C|) . Finally, in Section IVII 
we relate our experimental results with recent numerical 
simulations and draw some conclusions. 



II. THE EXPERIMENT 

A schematic diagram of the apparatus is presented 
in Fig. ^). The experimental set up consisted of 
a horizontal rectangular tray with dimensions (cc, y) = 
(180, 90)mm. It was connected to an electro-mechanical 
shaker so that a mixture of granular particles placed on 
the tray was vibrated longitudinally. Individual parti- 
cles or the granular layer were forced via stick-and-slip 
frictional contacts with the oscillating surface of the con- 
tainer. Different removable frames could be attached to 
the tray to make changes in size, shape and aspect ratio. 

The tray was mounted on a horizontal platform which 
was connected to a Ling LDS V409 electro-mechanical 
shaker. Its motion was constrained to be unidirectional 
by four lateral high precision linear bush bearings. The 
shaker was driven sinusoidally using a HP 33 120 A func- 
tion generator and the resulting dynamic displacement 
and acceleration of the shaking bed were monitored by a 
Linear Displacement Variable Transformer (LVDT) and 
a PCB quartz shear piezoelectric accelerometer. 

The main granular mixture consisted of phosphor- 
bronze precision spheres and poppy seeds and pho- 
tographs of representative particles are shown in Fig. ^) 
and d). The poppy seeds were non-spherical ("kidney" 
shaped) with an average diameter of LOGmm, polydisper- 
sity of 17% and a density of 0.2gcni~^. The phosphor- 
bronze spheres had a diameter of l.SOmrn, polydispersity 
of 3.0% and a density of 8.8gcm~^. In addition to being 
non-spherical the poppy seeds had a considerably larger 
surface roughness than the spheres, which resulted in a 
stronger frictional interaction with the surface of the os- 
cillating tray. This is evident from the series of ridges in 
the surface of the seeds, as shown in Fig. ^) 

We define the total filling fraction of the granular layer 
as, 

C[Nps, Npb) = — — — — = ipps + ippb, (1) 

xy 

where Nps and Npb are the numbers of poppy seeds (ps) 
and phosphor bronze spheres (pb) in the layer, Aps = 
(0.90 ± 0.15)mm^ and Apb — l.l&Tmm? are the two di- 
mensional projected areas of the respective individual 
particles and x and y are the longitudinal and trans- 
verse dimensions of the rectangular tray, ipps and Lppb 
are the individual filling fractions for each of the parti- 
cle species, respectively. Unless otherwise stated, C was 
varied by keeping the filling fraction of the spheres fixed 
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poppy seeds (pps- This enabled more controlled changes 
to be made in C but we also investigated changes in both 
ipps and Lppb, as will be discussed in Section IVBI 

In addition to the driving parameters and the total 
filling fraction, a geometrical dimensionless parameter, 
the aspect ratio, is defined as F = Ax/ Ay where Ax is 
the longitudinal length of the tray and Ay is the trans- 
verse width of the tray, perpendicularly to the direction 
of forcing. The values of Ax and Ay were changed using 
a variable frame which was positioned on all 4 side walls. 
Both Ax and Ay could be independently adjusted to the 
required value of F. 

All experimental runs were performed in an approx- 
imately monolayer regime. The larger heavy phosphor- 
bronze spheres were always in a monolayer but the lighter 
and fiatter poppy seeds could overlap. This degree of 
overlapping was due to both the polydispersity of the 
poppy seeds and the difference in size between the poppy 
seeds and the spheres (size ratio of g ^ 0.71). The mix- 
ture was deemed to be in the a monolayer regime if the 
extent of overlap of the poppy seeds was never such that 
the layer height exceeded the height corresponding to a 
diameter of the spheres. The failure of this criterion was 
readily noticed as smaller particles were observed to hop 
over domains of the larger particles. Hence, for the high- 
est values of the total filling fraction, C had values higher 
than those corresponding to maximum packing in two 
dimensions, which, for the case of monodisperse disks 
is 7r/-\/T2. This choice of performing the experiments 
in this approximately monolayer regime has two advan- 
tages. Firstly, particles are always in contact with the 
oscillatory surface of the tray, such that the forcing was 
provided homogeneously throughout the layer through 
frictional contacts. Secondly, these approximately two di- 
mensional experiments allowed the dynamics of the gran- 
ular layer to be fully visualized by imaging the system 
from above. 

All experimental runs were started using a homoge- 
neous mixture as initial conditions. This was achieved 
using the following procedure. Firstly, a particular fill- 
ing fraction of Nps poppy seeds were vibrated at large 
amplitude, A ~ ±5mm which created a homogeneous 
and isotropic layer. The phosphor-bronze spheres were 
then suspended above the layer, on a horizontal perfo- 
rated plate with (m x n) 2mm diameter holes arranged 
in a triangular lattice and held by a shutter on an inde- 
pendent superposed frame. The shutter was then opened 
and the Nph — m x n phosphor-bronze spheres fell onto 
the layer of poppy seeds, creating a near homogeneous 
mixture of the two types of particles. An example of 
such an initial configuration is presented in Fig. I^Ji). We 
found that if this procedure was not adopted, random 
initial clusters in a poorly prepared mixture could bias 
the results. 

The dynamics of the segregation process were visual- 
ized using a near homogeneous illumination of the tray 
and the behavior monitored from above using a CCD 
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high refiection coefficient compared to the poppy seeds 
and the spheres appeared as sharp bright regions so that 
direct observation of their motion was relatively straight- 
forward. The individual positions of the phosphor-bronze 
spheres were obtained in a central (73.1 x 59.5)mm^ vi- 
sualization window of the full granular layer in order to 
achieve the necessary resolution to obtain good estimates 
of the centers of the spheres. The spatial distribution of 
the positions of the phosphor-bronze spheres were ob- 
tained using image processing and particle tracking tech- 
niques which was then used to calculate a variety of mea- 
sures. No tracking was carried out on the poppy seeds. 

In summary, the control parameters of the experiment 
were the amplitude A and frequency / of the applied 
vibration, the total filling fraction of the mixture, C, and 
the aspect ratio of the container, F. 



III. THE SEGREGATION PROCESS 

In our previous work 0, we present experi- 

mental evidence for three qualitatively distinct phases of 
the binary mixture of poppy seeds and phosphor-bronze 
spheres. Each exists over a range of the filling fraction 
of the granular layer, C. We identify these phases as 
binary gas (unsegregated) , segregation liquid and segre- 
gation crystal |34| . 

The binary gas phase is found at low values of C and 
is essentially a collisional regime. In this phase there is 
enough free area and agitation, such that the particles 
move randomly around the layer and a mixed state per- 
sists. In the segregation liquid phase, at intermediate 
values of C, aggregation of the phosphor-bronze spheres 
occurs and mobile liquid-like clusters form. A representa- 
tive segregation pattern of the mixture in this segregation 
liquid phase is shown in Fig. ^p). The snapshot corre- 
sponds to a segregated state which self-organizes after 
the mixture is vibrated for a period of 3 minutes, with 
forcing parameters A — ±1.74mm and / = 12Hz and 
aspect ratio F = 2. The movement of the segregation do- 
mains is reminiscent of oil drops on water. The motion 
of the particles within the clusters is highly agitated and 
the collective motion is slow with merging and splitting 
of the domains. 

The transition from the binary gas to segregation liquid 
phase has the characteristics of a continuous phase 
transition with square-root dependence of the saturation 
levels as measured by the average stripe width which we 
treat as an order parameter. This critical phenomenon 
implies the existence of a transition point at Cc = 0.65, 
below which the layer remains mixed and above which 
segregation occurs. Moreover, critical slowing down of 
the segregation times scales is found near Cc- 

At high filling fractions, for C > 0.89 a second qualita- 
tive change of the structure and dynamics of the domains 
occurs as well defined striped patterns, perpendicular to 
the direction of forcing, self-organizes from the initial ho- 
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FIG. 2: Snapshots of the evolution of segregation domains in 
a mixture with C = 0.996, / = 12Hz, A = ±1.74mm (a) 
t = Qs - vibration of the granular layer was initiated from 
a homogeneous mixture, (b) t = 4.36s, (c) t = 9.36s, (d) 
t = 16.28s, (e) t = 23.2s, (f) t = 40.68s. The snapshots cor- 
responds to a central (73.1 x 59.5)mm'^ visualization window 
of the full tray. 



regation crystal. The monodispersity of the phosphor- 
bronze spheres means that the segregation domains in 
this phase consist of particles disposed in a hexagonally 
packed lattice whereas the polydisperse poppy seeds al- 
ways move around randomly. 

A time sequence which is typical of the behavior seen 
at at high filling fractions where crystalline stripes form, 
is presented in Fig. |21 for C = 0.996. Immediately after 
the vibration is applied, single large particles diffuse in 
a sea of the smaller ones, exploring different local con- 
figurations. When two large particles happen to come 
close together, the smaller particles cannot fit between 
them, and hence the pair is subjected to an asymmet- 
ric pressure that keeps it together. Subsequently, pairs 
may encounter others so that progressively larger clus- 
ters form. The unidirectionality of the driving induces 
an asymmetry in the segregated domains such that elon- 
gated domains of the larger phosphor-bronze spheres de- 
velop in a direction which is orthogonal to the direction of 
the drive. During this initial period, the rapid formation 
of clusters suggests an effective attractive force between 
the phosphor-bronze spheres, that leads to aggregation. 
Eventually, long domains form in the y-direction and well 
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self-organization into segregation domains occurs within 
timescales of tens of seconds. For longer timescales of 
the order of a few hours, these segregated domains pro- 
gressively coarsen with time, thereby merging to form 
increasingly robust stripes. The width of the domains 
follows a t^/** power-law with time and this scaling is in- 
dependent of the mixture used |32j | . The coarsening takes 
place most obviously for values of C in the range ~ 0.1- 
0.2 above Cc- An extensive parametric investigation of 
the coarsening behavior as a function of C has yet to be 
carried out. 

It is important to stress once more that, depending 
on the control parameters, the mixture does not always 
evolve into a robust striped pattern. This was only for 
filling fractions with C > 0.89. At lower filling fractions, 
in particular closer, but above, to the segregation tran- 
sition point Cc the segregated domains are increasingly 
mobile and blob-like, as shown in the representative liq- 
uid state of Fig. 



IV. VORONOI MEASURES 

We now focus on the description of the segregation in 
our granular mixture using a number of microscopic mea- 
sures. By microscopic, we mean that both structural and 
dynamical quantities are analyzed using the positions of 
the centers of the individual phosphor-bronze spheres ob- 
tained from the particle tracking analysis. This is by way 
of contrast with the macroscopic average width of the seg- 
regated domains and respective fluctuations used in our 
previous work [SJ, yj, |33 . 

The construction of Voronoi cells through tessellation 
(also know as Wigner-Seitz cells) is a standard tool for 
the study of spatial configurations of particle ensembles 
which is widely used in condensed matter physics '36*1 and 
is outlined as follows. Consider a set P of coplanar par- 
ticles with their centers located at Ci{x, y) for i = 1 — > N. 
For each particle i, Voronoi tessellation yields a polygo- 
nal cell that encloses a region inside which any point is 
closer to the center Ci of the i"* particle than any other 
in the set P. We have used the voronoi (x,y) routine in 
the package MATLAB 7.0 to implement this procedure. 
It is straightforward to extract a measure of the local area 
density associated with each phosphor-bronze sphere for 
this geometrical construction. It is also possible to ob- 
tain another useful measure of the angle between nearest 
neighbors. These two quantities are introduced in the 
following two Sections. IIV Al and llV Bl respectively. 



A. Local Voronoi density 

The first quantity considered is a measure of the lo- 
cal area density associated with each phosphor-bronze 
sphere . Following a standard procedure |33| , the lo- 
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video frame can be defined as the ratio, 

-^sphere 



Pv 



^cell 



(2) 



where Asphere = 7r((i/2)^ is the two-dimensional pro- 
jected area of the imaged spheres with diameter d and 
Al^ii is the area of its Voronoi polygon. 

In Fig. 121 we present examples of typical Voronoi con- 
figurations, at four different values of C, constructed us- 
ing the positions of the phosphor-bronze spheres. In the 
binary gas regime at C = 0.495 no segregation occurs 
and the network of Voronoi polygons appear random, as 
shown in Fig. ISfa). By way of contrast, at C = 1.007 
where definite segregation develops, structure appears in 
the domains as can be seen in Fig. Old). Two snapshots 
of the network of Voronoi cells for intermediate values 
of C are presented in Fig. ^) and (c). Note that the 
particles at the edges of the segregation clusters have an 
associated area density significantly lower than those in 
the bulk of the domains. 

We first discuss the time evolution of the local Voronoi 
area density from the initial mixed state. A time window 
of At = 4:Sec, which corresponds to 100 video frames 
(i.e. ~ 48 drive cycles) has been used to obtain dynamic 
averages for the area density of individual spheres, p*, 
as, 



Pv 



7iAt 



(pl) 



(3) 



where the brackets (.)„ denote averaging over all the par- 
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FIG. 3: Voronoi diagrams obtained from the positions of the 
phosphor-bronze spheres, for binary mixtures with various 
filling fraction values: (a) C = 0.495, (b) C = 0.623, (c) 
C = 0.687 and (d) C = 0.751. The frames correspond to 
configurations obtained 40sec after vibrating an initially ho- 
mogeneous mixture. 




100 150 
Time (s) 

FIG. 4: (a) Time evolution of mean local Voronoi area density, 
having started from initially homogeneous mixtures. (♦) C = 
0.495, (•) C = 0.687, (■) C = 1.049. 



tides, z, found within the rt*'^ time window, rtAr < i„ < 
(n -I- l)Ar with n G [0, 1,2,3, ..,59]. A time dependent 
Probability Distribution Function for the local Voronoi 
area density, PDF{^, t„), was obtained by constructing 
normalized histograms of as a function of the discre- 
tised time, i„. Each PDF(j)^,tn), for any time window 
n, typically contained statistical ensembles with 35,000 
to 40,000 particles. 

In Fig. 21 we plot the mean value of the distribution 
t„), Dy{tn), for three values of the filhng frac- 
tion. Let us first focus on the behavior at early times, i.e. 
within the period t < 100s. At C = 0.495, £>„(*„) is fiat 
since the layer remains mixed. At intermediate values of 
C, where C = 0.687 is a typical example, Dy{tn) exhibits 
a slow increase up to a value of Dy{tn) ~ 0.3 as segre- 
gation clusters form. At high values of filling fraction, of 
which C = 1.049 is representative, there is a rapid ini- 
tial evolution, since increasingly dense clusters form, up 
to a value of Dy{tn) ^ 0.525 after which the mean area 
density levels off. This behavior is consistent with the 
observation of saturation in the macro scop ic mean stripe 
width presented in our previous work |33l |3J, |35| . 

At later times, i.e. t > 100s, the mean Voronoi area 
density for the two cases of C = 0.687 and C = 1.049 
only shows small deviations from the level off value of 
Dv{tn) and the system has reached a segregated state 
with a characteristic Voronoi density. The small dips in 
Dv{in) correspond to the long term splitting and merging 
of stripes mentioned in Section UllI An error estimate of 
the closeness to a steady state can be obtained from the 
ratio of the standard deviation to the mean of Dy{177 < 
t < 233sec), which was below 6.5% for all values of filling 
fraction considered. 

The dynamical behavior of the segregation process can 
now be analysed as a function of the filling fraction using 
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of the formation of domains. During this regime of fast 
initial segregation growth, Dy exhibits an approximately 
linear behavior of the form, 



Dy{C,t) = 0) + 6(C). < 



(4) 



where Dy(t = 0) is the mean area density of the initial 
homogeneous mixture and 6(C) is the corresponding rate 
of segregation. A value of Dy{t = 0) = 0.158 ±0.004 pro- 
vides a good fit to all experimental runs since the layer is 
consistently started from an homogeneous mixture. The 
inverse of the rate of segregation yields a measure of the 
segregation timescale, toiC) = 1/6(C). The quantity 
tniC) is plotted in Fig. El The segregation timescale 
rapidly increases as Cc is approached from above, indi- 
cating a slowing down of the dynamics which is consistent 
with the critical slowingdown found previously from the 
macroscopic measures |33l Issj ] . 

In 34] we showed that the behavior of the PDF{py) 
distributions obtained in the steady state regime (i.e. af- 
ter the initial segregation growth) is useful in the charac- 
terization of the segregation in a mixture. At low values 
of C, the PDFs are peaked at small pv As C is increased 
a qualitative change in the shape of the PDFs is seen and 
at C '--^ 0.65 they flatten out indicating that there is a 
greater probability of finding particles with an area den- 
sity across the entire range. As C is increased further, a 
new peak develops at high area densities corresponding 
to particles within the segregation clusters. This peak at 
high py becomes increasingly sharper for high C, with a 
drop at values of p„ — 0.9, which is consistent with max- 
imum packing in 2D of t:/^/V1 for a perfect hexagonal 
arrangement of disks. 

The value of the median of the PDF[py) distributions, 
p^°'^ , measures the characteristic Voronoi density of the 
spheres in the mixture and is used in Section to aid 
mapping out the phase diagram of the granular mixture. 
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FIG. 5: Segregation timescale, to ~ l/b{C) of the phosphor- 
bronze spheres as a function of filling fraction. The solid 
vertical line is positioned at the critical point for segregation, 
Cc obtained in 1341. 
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FIG. 6; Probability distribution functions for the angle be- 
tween nearest neighbors, PDF [5): for (•) C = 0.516, (♦) 
C = 0.687, (A) C = 0.815, (■) C = 1.028. Inset: Schematic 
diagram for the definition of angles between nearest neigh- 
bours, 5. A, B and d are the positional coordinates of three 
neighbouring particles. The angle defined by the three parti- 
cles, about Ci, is 5i = ai + Pi. The solid lines are a section 
of the Voronoi polygon, belonging to the particle d, which 
is defined by the vertices Vi, V2 and V3. 



B. Angle between nearest neighbors 

Another quantity that can be calculated from the 
Voronoi tessellation procedure is the angular distribution 
between nearest neighbours of each sphere. The configu- 
ration of a section of a schematic Voronoi polygon for a 
particle with coordinates C; and two of its nearest neigh- 
bours, A and B, is given in the inset of Fig. |S| We define 
Si to be the angle between nearest neighbours set by A, 
B and C^. This can be obtained from the Voronoi tessel- 
lation procedure which yields the position of the coordi- 
nates of the vertices of the Voronoi polygon Vi , V2 and 
V3. Hence, it follows from the geometric construction 
around the particle that, 



5, = + = 180- (0,+^,) 



where. 



and 



cos 



(Vi -V2).(C,- V2) 
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Representative curves of PDF{5) are presented in Fig. 
inifor four different values of C. For example, for small C 
the distribution is broad which indicates that the spheres 
do not exhibit a preferred orientation as expected for a 
mixture. However, as the filling fraction is increased the 
main peak of the distribution shifts towards 60° and the 
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The location, 5"*°^, and width, w, of the peak of the 
PDF{5) distributions is therefore a measure of the lo- 
cal orientational structure nearby each phosphor-bronze 
sphere. Both of these quantities are used in Section 
to aid mapping out the phase diagram of the granular 
mixture. 



V. EXPLORATION OF PARAMETER SPACE 
A. Aspect ratio 

Since C is a dimensionless ratio between the total pro- 
jected area occupied by all the particles and the area 
of the tray, it is now of interest to explore the effect of 
varying the total area and aspect ratio of the tray. A 
frame with movable walls in both x and y was used to 
accomplish this. Its side walls were mounted in preci- 
sion machined parallel grooves and the respective dimen- 
sions adjusted using a Vernier scale which was accurate 
to within 0.05mm. The range of the aspect ratio explored 
was 0.22 < r < 9.00. The filling fraction of the layer was 
fixed at C = 0.900 with t/jpf, = 0.174 phosphor-bronze 
spheres. Since the total area of tray varied while chang- 
ing r, both Nps and Npi, had to be changed accordingly in 
order to keep C and ip constant. The forcing parameters, 
as before, were set at / = 12Hz and A = ±.1.1 Aram. 

Example snapshots of the granular layer taken after 
2min of vibration of an initially homogeneous mixture 
are presented in Fig. [7{e-f), for various values of F. Note 
that the photographs correspond to a central imaging 
window of the full tray (73.1 x 59.5mm^) and hence in 
all experimental frames of Fig. |7| only a portion of the 
granular layer is shown. The vertical boundaries of the 
tray are, therefore, only observable in Fig. [7| (a) and (b) 
and the horizontal boundaries in Fig. [7| (e) and (f ) . 

For F <C 1, of which the frame in Fig. [7{a) is a rep- 
resentative example, a single stripe formed perpendicu- 
lar to the direction of forcing. For larger values of F, a 
number of similar well defined vertical stripes rapidly de- 
veloped in the same way of the particular case of F = 2, 
which was studied in detail in Sections IIIII IIV Bl More- 
over, the number of stripes in the system increased as F 
was increased. For example, increasing the tray's aspect 
ratio from F = 0.222 (Fig. |7|i) to F = 0.667 (Fig. {7}p) 
resulted on a increase from one to two stripes. 

Even if larger number of stripes could be attained by 
varying F, in all cases the stripes had a similar form i.e. 
they had approximately the same longitudinal width even 
if their transverse height (set by Ay of the tray) decreased 
with F. The characteristic local Voronoi area density of 
the spheres, p™"^ was measured as a function of F. As 
the aspect ratio of the tray was changed, the local area 
density of the phosphor-bronze spheres remained approx- 
imately constant (to within 6%) with an average value of 
(p™°^)r = 0.696 ± 0.040. This suggests that the segre- 
gation behavior appears to be independent of both the 
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FIG. 7: Snapshots of the segregated phase for various values 
of F. The images correspond to the layer configuration, 2min 
after vibrating an initially homogeneous mixture with the fol- 
lowing parameters: C — 0.900, (fipt = 0.174, / — 12Hz and 
A = ±1.74mm. (a) F = 0.222 (b) F = 0.667 (c) F = 1.556 
(d) F = 2.571 (e) F = 4.500 (f) F = 9.000. The end-wall 
boundaries can be seen as dark vertical bands in (a) and (b) 
and the sidewalls as horizontal bands in (e) and (f). 



B. Filling fraction parameter space 

The majority of the experiments were performed 
by keeping the filling fraction of the phosphor-bronze 
spheres fixed at tpph — 0.174, and incrementally increas- 
ing the filling fraction of the poppy seeds, tpps. Hence it 
is possible that the transition sequences might be caused 
by reaching a critical fill ratio of the domain. In or- 
der to make sure that this was not the case, a series 
of experiments were performed where the dependence of 
the segregation transition on the relative composition of 
the binary mixture was investigated, while keeping the C 
constant. The amounts of spheres and poppy seeds, set 
by ippb and (pps, were varied accordingly, with the con- 
straint of keeping the total filling fraction of the layer 
fixed at C = 0.79 ± 0.06. The mean longitudinal width, 
a, of the domains of phosphor-bronze spheres has been 
used as the order parameter to measure the state of the 
system. Details concerning this measure can be found in 
[s^,!!^. A plot of a versus ippb is given in Fig. |S1 where it 
can be seen that for low values of ippb, the mean domain 
width is cr = 3.23 ± 0.04. This corresponds to approxi- 
mately two sphere diameters and is consistent with the 

phfinrp nr-piirrpripA nf npitrlilinrin0' nnrtiplps PArp.rp.crpr] nv(=>r 
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FIG. 8: Mean stripe width of spheres as a function of tpph- 
Total filling fraction was kept constant at C = 0.79 ± 0.06. 



the layer. Therefore, these are not segregated domains 
and the layer is in a mixed state. As ^pph is increased 
past (pp^=0.055, cr increases continuously. We emphasize 
again that tppg has to be changed accordingly in order to 
keep C constant. This result suggests that, for segrega- 
tion to occur, not only has the total filling fraction to be 
lar ge e nough (i.e. above a critical value Cc as reported 
in [231) but also the amount of spheres in the layer is 
required to be above a threshold value ip^^^. This result 
emphasizes the fact that the total filling fraction is a two 
dimensional parameter C{ippb, ipps)- 

Therefore, the dependence of the segregation behavior 
on C was further explored in the filling fraction param- 
eter space, {ifps, fpb) by increasing the number of poppy 
seeds and hence scanning along (ppg for four values of 
<fpt = {0.054,0.083,0.118,0.174}. 

The measure used to characterize the state was Voronoi 
area density, p™"^, and the results are presented in Fig. 
IHl^a-d). All four data sets exhibit a clear transition be- 
tween a mixed state (low area density) and segregated 
(high area density) state. A particularly interesting fea- 
ture is that the transition becomes sharper in the mix- 
tures with lower values of ippb and yet the position of 
the segregation critical point remains fixed. The mean 
value for the estimate of the mixed state at the left hand 
edge of each graph is proportionally lower as (fpb is de- 
creased. This is expected since the overall area density of 
the phosphor-bronze spheres in the mixture is decreased, 
i.e. the total number of spheres in the tray is smaller for 
lower ippf). However, the branches corresponding to the 
segregated phases do not show any significant change, for 
different ifpb- Estimates of the local densities of ~ 0.8 
were attained in all four cases, at the highest values of 
C. 

The location of the maximum of the PDF{5) distribu- 
tions presented in Section FlV Bl yields the characteristic 
angle between nearest neighbors, ,5™°^. In Fig. IHIe-h) we 
present the results of (5™°^ for the four data sets with dif- 



becomes apparent in this measure for lower values of ippb- 
Take, for example, the data in Fig. |^) for ippb = 0.054. 
At the lowest values of C, J™"^ tends to 45°. However, 
as C was increased, J™"^ follows a non-monotonic varia- 
tion, as the mixture goes through the segregation transi- 
tion. The characteristic angle first decreases to 38° only 
to then sharply increase past the segregation transition 
point. The nature of this non-monotonic dependence of 
gmax (J jg unclear, even though it appears to be robust 
and persists in the other three datasets. At the high- 
est values of C the phosphor-bronze spheres self-organize 
into a nearly hexagonal pack, for which one would ex- 
pect SJ^J!f = 60°. This is consistent with the value of 
gmax ^ (-57 2 ± 2.0)° attained in the experiment at high 
C. 

The second quantity that we extract from the PDF{5) 
distributions is the width at half maximum, w, of the 
peak which is a measure of the fluctuations in S and is 
plotted in Fig. IHIi-l), for the four values of (ppb- In all 
curves, three distinct regimes are evident where each cor- 
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FIG. 9; Voronoi measures for the filling fraction phase space. 
First column: (a-d) Characteristic Voronoi area density. Sec- 
ond Column: (e-h) Characteristic angle between nearest 
neighbours. Third Column: (i-1) Width of of the peak of 
the PDF(5) distribution. The x-axis of all plots is the To- 
tal filling fraction of the layer, C with ippt set to (a,e,i) 
ifpt = 0.054, (b,f,j) ippb = 0.083, (c,g,k) ifpt = 0.118 and 
(d,h,l) ifipi = 0.174. 
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FIG. 10; Phase diagram for the {ippb,(pps) parameter space, 
(o and A) binary gas. (□ and *) segregation liquid. (+) 
segregation crystal. The points along the oblique solid line, 
1/2, correspond to the data presented in Fig. ISlfor which the 
total filling fraction was kept constant at C = 0.79 ±0.06 and 
the relative amounts of particles were changed accordingly 
(through ippt and ^ps). The dashed lines, LI and L3 are the 
experimental phase boundaries between the three phases of 
the mixture. 



responds to one of the phases of the binary mixture. As 
an example, consider the curve for ipph = 0.118 given in 
Fig. Eli). At low C the distribution is broad, indicating 
that a wide range of angles between nearest neighbors 
are possible. This implies that a disordered binary gas 
phase exists and there is no segregation. At Cc ^ 0.7, the 
point at which the segregation transition occurs, there is 
a quantitative change of behavior and w rapidly drops 
with increasing C. This region corresponds to the seg- 
regation liquid phase. At Cm ~ 0.9, w levels off at low 
values and PDF {5) is now sharply centered &t 5 — 60°, 
implying that the transition from the segregation liquid 
into the segregation crystal has taken place. This three 
step behavior in w can be seen in all four datasets con- 
firming that the scenario of three segregation phases as a 
function of C is robust over a range of the filling fraction 
parameter space. 

In Fig. El the four data sets discussed thus far are 
combined in a {(Pps, Vpfc) phase diagram where the phase 
boundaries Cc and Cm were calculated from the points 
demarking the three regimes of w discussed above. The 
aim is to determine the locus of existence of the three 
segregation phases. The binary gas is located in the left 
hand side of the diagram in the regions of low (pps . The 
segregation liquid exists in the central region of the dia- 
gram. The segregation crystal is observed at large values 
of ipps , on the right hand side of the phase diagram. Note 
also that the phase boundaries LI and L3 (dashed lines 
in Fig. I10|) are approximately parallel for the data sets 
with cppb = {0.083,0.118,0.174}, to within - 9%. This 

inrlipfi.tps that in this rpcrinn nf th(=> nhnsA rlincrrpi.Tn thp. 



total filling fraction is the primary parameter in deter- 
mining the behavior, rather than the absolute amounts 
of each of the particle types. However, larger deviations 
occur for the dataset with ippi, — {0.054}. The devia- 
tions in this region of the phase diagram with low ippb is 
confirmed by the dataset discussed above where the total 
filling fraction was kept constant at C = 0.79 with (pph 
and ifps varied accordingly (data along the oblique solid 
line in Fig. . 

We stress that strong deviations from the results pre- 
sented here are to be expected in the limiting cases as 
the mixture approaches the single component regimes: 
i) 'Ppb and ii) ipps — > 0. Two sub-limits of partic- 
ular interest are those for the dense cases: iii) {(ppb — > 
y/l2/Tr, ipps 0) where crystallization of the phosphor- 
bronze spheres prevails and {ippb 0, (pps ^1) where 
the dense liquid- like nature of the poppy seeds dominates 
and segregation is suppressed. This is consistent with the 
polydispersity and shape of the poppy seeds preventing 
crystallization. 



C. Forcing parameter space 

In the results discussed thus far, the drive parameters 
have been kept fixed at / = 12LIz and A = ±1.74mr7i. 
We now report the results of an investigation of the (/, A) 
parameter space. The aspect ratio and total filling frac- 
tion were fixed at F = 2 and C = 0.900 with (ppb = 0.174. 
All experimental runs in this section, as before, were ini- 
tiated from a homogeneously mixed configuration. 

With {f,A) = (12 Hz, ±1.7 4mm) the above conditions 
produced a dense segregation 'liquid' from the initial mix- 
ture. However, all three phases can be obtained at ap- 
propriate locations in the (/, A) space. Examples of each 
of these phases can be seen in the experimental frames 
presented in Fig. Ilir a-c). An example of a binary gas is 
shown in Fig. [nja) for (/, A) = (28ffz, il.lSmm). The 
forcing is sufhciently large to induce apparently random 
motion in the phosphor-bronze spheres across the layer 
and no segregation occurs. This example of a binary gas 
differs from the cases considered previously for C < Cc in 
the sense that, at this high filling fractions (C — 0.900), 
the particles in the layer have persistent contacts with 
their neighbors and there is only a small amount of free 
area in the tray. Therefore this is a highly agitated but 
non-coUisional state. At {f,A) = (IGTJz, ±1.50mm), 
segregation occurred and the segregation liquid shown 
in Fig. Illf b') was observed. A representative exam- 
ple of a segregation crystal is shown in Fig. Illf c) for 
{f,A) — (12i?z, ±0.58mm), where spheres within the 
domains were disposed in hexagonally packed configu- 
rations. 

The corresponding phase diagram for the (/, A) pa- 
rameter space is presented in Fig. IllH ). Two clear ex- 
amples of the segregation crystal phase were found at 
parameter values located in the lower right hand corner 

nf th(=> rlifi.crrnrn i p fnr Inw amnlitnrlps p.nrl frpmipncips 
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At intermediate values of the forcing, a segregation liquid 
phase was observed. At relatively large values of / and 
A segregation did not occur and a binary gas phase was 
dominant. The nature of the phase boundary between 
the segregation liquid and the binary gas phases (solid 
line in Fig. Illl i) is discussed below. 

In order to obtain quantitative estimates of the behav- 
ior, the PDF{pv) of the distribution of the local Voronoi 
area density was calculated for each point in (/, A) space. 
The characteristic local Voronoi area density, p™""^, was 
then obtained from the location of the maximum of the 
distribution, as before. 

The dimensionless maximum acceleration of the tray 
is used to parameterize the forcing. It is given by , 



7 = 47r' 



,Ap 



(8) 



where g is acceleration due to gravity. The parameter 
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FIG. 12; Plot of the characteristic local Voronoi area density 
of individual spheres versus the non-dimensional maximum 
acceleration, 7 for a range of values of (/, A). The transition is 
between a mixed (horizontal) and segregated (sloped) states. 
The solid and dashed lines are the best least squares fit in the 
segregated and mixed regimes, respectively. The point with 
the lowest 7 corresponds to a layer for which the forcing was 
insufficient to agitate the granular mixture. 
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FIG. 11: Experimental frames of the binary mix- 
ture at various values of / and A. (a) Binary 
gas, {f,A) — (28-ff2, ±1.18mm). (b) Segregation liq- 
uid, (/, A) = (16-ff2:, ±1.50mm). (c) Segregation crystal, 
(/, A) = (12ii'z,±0.58mm). F = 2 (Ay = 180mm) and 
C = 0.900 with ifpb = 0.174. (d) Phase diagram for (/,A) 
parameter space showing regions of existence for the three 
phases, (o) Binary gas. (□) Segregation liquid. (-I-) Segre- 
gation crystal. The points A, B and C, along the parameter 
path r, correspond to the experimental frames shown in (a), 
(b) and (c), respectively. The phase boundary between the 
segregation liquid and the binary gas phases is represented by 
the soUd black curve and is given by Eqn. @ with 7c = 2.95. 



7 is commonly used in vertically vibrated granular sys- 
tems |38| in which gravity plays a dominant role and the 
granular layer requires 7 > 1 to leave the vibrating base. 
In our horizontal set up, gravity enters the problem in- 
directly through the frictional forces acted on individual 
particles, F = pmg where m is the particle mass and /i 
its friction coefficient. In the ideal scenario of no rolling, 
the value of the non-dimensional acceleration at which 
relative motion between the particle and the oscillatory 
tray appears would occur at 7 = /i. 

In Fig. ^] the characteristic local Voronoi area den- 
sity is plotted as a function of 7. All data points from 
the (/, A) parameter space study collapse on to a sin- 
gle curve, where two clearly distinct regimes can be ob- 
served. This indicates that 7 is, indeed, an appropriate 
parameter to describe the forcing of the granular mixture 
through oscillatory motion of the tray. 

At high accelerations, for 7 > 3, the characteristic 
Voronoi area density remains approximately constant 
with increasing 7. In this regime, the layer was in 
the highly agitated binary gas phase discussed above 
and no segregation occurred. As 7 is decreased below 
7 '--^ 3 the final state achieved becomes increasingly more 
dense in an approximately linear fashion. The solid and 
dashed lines in Fig 1 121 were obtained from the best least 
squares fits of the form p™°^(7) = m7 -I- n, in the region 
0.34 < 7 < 2.67 (segregation phase) and /9™°^(7) — p 
in the region 3.34 < 7 < 5.37 (binary gas phase), re- 
spectively. The intercept of the two lines yields the lo- 
cation of the transition point which was measured to be 
^_ = 9 Q.c; -I- n 1 fi 
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The relationship, 




provides a good fit to the phase boundary, f'{A'), be- 
tween the binary gas and segregation hquid phases as 
shown by the sohd curve in Fig. ITTH). This also sug- 
gests that 7 is the appropriate parameter to describe the 
forcing of the tray. 

VI. CONCLUSION 

We have carried out a detailed experimental investiga- 
tion of segregation behavior in a monolayer of two types 
of particles contained in a horizontal oscillating tray. The 
particles are set into motion via stick slip interaction 
with the base and this effectively randomizes their mo- 
tion. An initially homogeneous mixed layer gives rise to 
spontaneous robust patterns when set into motion. The 
patterns have the form of clusters of one particle type 
and they persist over a wide range of the experimental 
conditions. The two principal control parameters of the 
system were found to be the layer filling fraction and the 
dimensionless acceleration of the driving. Moreover, we 
demonstrated that the essential mechanisms of the segre- 
gation process appeared to be independent of the aspect 
ratio of the tray. 

Three qualitatively distinct phases were identified: a 
disordered binary gas for low values of C, segregation 
liquid with mobile domains of one of the particle types 
at intermediate values C and segregation crystal where 
the domains are stripes oriented in a direction orthogonal 
to that of the driving at high values of C. Moreover, we 
have reported a novel transition between the segregated 
and mixed phases as the dimensionless acceleration of the 
driving was increased. 

Recently, there has been a number of numerical stud- 
ies where systems analogous to that of these experiments 



have been simulated [s^ 0, 0, 0. The segregation 
behavior of these numerical systems is found to be in 
partial qualitativ e ag reement with our experiments. Cia- 
marra et. al. [4lll42l | have suggested that the segregation 
process is a result of dynamical instability which resem- 
bles the classical Kelvin-Helmholtz instability observed 
at a fluid interface. Similar stripe formation has also 
been observed in a continuum model of a binary fluid in 
which the components are differentially forced js^ . 
These studies suggest that the principal mechanism un- 
derlying segregation in this class of systems is the dif- 
ferential driving between each of the particle species. In 
our experiments, this could relate to the different surface 
properties of the phosphor-bronze spheres and the poppy 
seeds which induce different frictional interactions with 
the surface of the oscillatory tray. 

However, even though the segregation patterns of the 
simulations are qualitatively similar to those observed 
in our experiments, none of the numerical models are 
as yet able to reproduce the critical behavior. In par- 
ticular, there is no evidence for the associated critical 
slowing down of the dynamics near the transition region. 
Moreover, our parametric study of the driving parame- 
ters and the additional transition we have uncovered as 
a function of the dimensionless acceleration of the tray 
demonstrates that the forcing also plays a crucial role in 
the phase behavior of the granular mixture. This point 
has so far been overlooked in the numerical studies. We 
believe that this extensive investigation of the parameter 
space of our experimental system will be crucial for fur- 
ther testing and refinement of numerical models of the 
type mentioned above. 
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